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Abstract. The question of complete integrability of evolution equations associated 
to n X n first order isospectral operators is investigated using the inverse scattering 
method. It is shown that for n > 2, e.g. for the three-wave interaction, additional 
(nonlinear) pointwise flows are necessary for the assertion of complete integrability. 
Their existence is demonstrated by constructing action-angle variables. This con- 
struction depends on the analysis of a natural 2-form and symplectic foliation for the 
(~| . groups GL{n) and SU(n). 
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More generally, k independent commuting Hamiltonian flows in a 2A^-dimensional 
manifold are said to be a completely integrable family if there are N—k independent 
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^ ■ 1. Introduction 

(N 

CN I A classical Hamiltonian flow with 2N degrees of freedom is said to be completely 

t^^ , integrable if it has A^ independent integrals of the motion which are in involution. 

-(— > 
c^ 
S I integrals of the motions which are in involution, or equivalently the N — k flows 

^ ! may be enlarged to a set of A^ independent commuting flows. By a theorem of 

S^ . Jacobi and Liouville, there then exist (at least locally in phase space) a new set of 

c^ ■ canonical variables, called action-angle variables, in which the flows are particularly 

simple; set [A] for a precise global version due to Arnold. 

In recent years a number of nonlinear evolution equations, beginning with the 
KdV equation, have been shown to have Hamiltonian form on appropriate inflnite- 
dimensional manifolds and to have an inflnite family of integrals of the motion 
which are in involution. Such equations are commonly referred to as "completely 
integrable" , although it no longer makes sense to count half the number of dimen- 
sions. Nevertheless the inverse scattering method makes it possible to give a precise 
form to the question of complete integrability and, indeed, to reduce it to a question 
in a flnite dimensional space. 

Such results are known for the KdV hierarchy and the nonlinear Schrodinger 
hierarchy; [Ga] , [ZF] . The scattering map both linearizes and decouples these flows, 
and action angle variables have been obtained [ZF], [ZM]. In both these cases the 
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pointwise dimension of the scattering data is 2, and no new pointwise flows are 
needed to get to half the dimension. Manakov [Ma] obtained action-angle variables 
for the 3-wave interaction equation. 

The inverse scattering method is based on analysis of an associated linear spectral 
problem, and the associated hierarchy of flows are isospectral for the linear operator. 
For KdV the linear operator is the 1-dimensional Schrodinger operator and for NLS 
it is the 2 X 2 AKNS-ZS operator 

(1.1) ^-zJ-q{x); 

here z is the spectral parameter, and J is a given constant matrix. In this paper 
we consider the n x n version of (1.1), under the assumption that J is semisimple 
with distinct eigenvalues and that q takes values in the range of ad J and vanishes 
rapidly at oo. The three-wave interaction is an example of an associated isospectral 
flow, when n = 3; [ZM], [Ma], [Ka]. The scattering and inverse scattering theory 
for (1.1) has been considered by a number of authors: see [ZSl], [AKNS] for n = 2 
and [Ma], [Ka], [Sh], [Ne], [Ge], [BY], [BCl], [Ca]. 

Each traceless matrix fj, which commutes with J generates a hierarchy of isospec- 
tral flows of (1.1). On the scattering side these form an [n — l)-parameter family of 
commuting pointwise fows. The pointwise dimension of the space of scattering data 
for (1.1) is n^ — n. We show that there is an appropriate Hamiltonian structure 
on this space of pointwise data, and that the (n — l)-parameter family is com- 
pletely integrable in the classical sense: it is part of an (n^ — n)/2-parameter family 
of commuting Hamiltonian flows. Both the existence of an appropriate pointwise 
Hamiltonian structure and complete integrability follow from a construction of Dar- 
boux coordinates (coordinates which diagonalize the 2-form) which are action-angle 
variables for the flows of the hierarchies. It should be noted that the additional 
commuting point-wise flows needed when n is greater than 2 are not linear on 
scattering data. 

The results just described are obtained in the category of complex manifolds and 
Hamiltonian structures. We are also interested in the real case. The three-wave 
interaction, for example, is associated to the operator (1.1) with J+J* = and the 
constraint q + q* = 0; on the scattering side the appropriate group is SU{3) rather 
than SL{3, C). We show that one can flnd real Darboux coordinates for scattering 
data to provide action-angle variables for the 3-wave interaction and the other 
flows of the hierarchies. The canonical transformation to action-angle variables is 
not algebraic in this case: it requires the Liouville method and elliptic functions. 
Manakov [Ma] used a different method to obtain action-angle variables for the 3- 
wave interaction which have a simple form but which are nonlocal functions of the 

— j-«;„„ „f j-u„ „„„j-j- — ; — , »„„j-„;,.. „ „f „ — .j-; — o. j-U„ „„„„„;„j- — 1 a „ . „i„ i„„„i 



in s. 

Our analysis of the Hamiltonian structures (symplectic form, Poisson brackets) 
leads to a natural closed 2-forni of rank n^ —n on GL{n), and a natural symplectic 
foliation of GL{n). The reduction J + J* = 0, q' + Q'* = leads to consideration 
of SU{n) in place of GL{n). The induced Poisson bracket is not a Poisson-Lie 
structure [Dr], since it is not degenerate at the identity element. However it was 
pointed out to us by Lu [Lu] that our structure is the translate by a Weyl element 
of a Poisson-Lie structure which is the classical limit of a quantum group structure 
described by Drinfeld [Dr]. 

The plan of the paper is the following. In section 2 we review the Hamiltonian 
structure and hierarchy of flows associated to the operator (1.1). The scattering 
theory for the case J+J* = is reviewed in section 3. We then compute the Poisson 
bracket for scattering data and state the main results on existence of Darboux 
coordinates and complete integrability. In section 4 we introduce and analyze the 
2-form on GL{n) and obtain Darboux coordinates. A symplectic foliation of GL{n) 
is introduced in section 4, and we calculate the associated Poisson bracket and the 
Hamiltonians for a family of linear flows. 

The algebraic results of sections 4 and 5 are used in section 6 to prove the results 
on Darboux coordinates and complete integrability for scattering data which were 
stated in section 3. The case of SU{?>) is taken up in section 7; complete integrability 
of the three- wave interaction is a consequence. In section 8 we show that the results 
stated in section 2 remain valid valid without the restriction J -|- J* = 0. 

We are grateful for discussions with V. Zurkowski, J. H. Lu, G. Moore, and G. 
Zuckerman. B. Konopelchenko drew our attention to the work of Manakov. 

2. Symplectic structure of Hamiltonian hierarchies 

We consider Hamiltonian hierarchies of flows associated to the first order differ- 
ential operator 

(2.1) ±^-zJ-q{x), zee, 

where J is a constant nxn semisimple matrix; q{x) is an n x n matrix whose entries 
qjk belong to the Schwartz class 5(R); and, for each x, q{x) lies in the range of 
adJ. We denote by P the linear space of all such q; thus P = 5(R; a(iJ(M^)), 
where M^ = M„(C) is the space of n x n matrices, with the Schwartz topology. 
We use the following inner product on P: 

(2.2) (q^p) = / tr\q(x)p(x)]dx. 



Since P is a linear space we may identify it with its tangent space. We denote 
tangent vectors (at a given point q) by q. Associated with (2.1) is the closed 2-form 

(2.3) np = \ f tr[5q{x) A [adJ]-^5q{x)]dx 

where 6q{x) denotes the linear functional taking q to q{x) and [adJ]~'^ maps to the 
range of adJ , on which adJ is injective. Thus 

(2.3') ^p{qiA2) = - I tr[qi{x)[adJ]~^q2{x)-q2{x)[adJ]~^qi{x)]dx 

^ Jr 

Since the inner product is non-singular, Op is symplectic. Note that when J+J* = 

0, and we restrict to the set {q & P : q + q* =0}, then the form Op is real. 

We shall work with the case in which J is diagonal with distinct eigenvalues: 

J = diag(zAi, ■iA2, . . . iXn) ■ In this case 

(2.4) Op = / y2~r\ :^-^Sqjk{x) A6qkj{x)dx. 

A Poisson bracket is associated to the symplectic form Op in the standard way. 
If F is a functional on P which is Frechet differentiable, and g is a tangent vector, 
we write 

(2.5) [qF]{q) = ^ |^^^ F{q + eq) = (^(?),^), 

i.e. we identify dF/dq with the gradient of F. The Hamiltonian vector field associ- 
ated to F , denoted Hp, is then defined by 

(2.6) np{HF,q) = -qF 

This definition and (2.3') imply that Hp = [J,dF/dq]. The Poisson bracket [Ne] is 
then given by 



(2.7) 



{F,G}p = HfG = -Op(ifG,^F) = J tr ([J, ^-f]^-f] dx. 



There is an (n — l)-parameter family of hierarchies of commuting Hamiltonian 
fiows in P, defined as follows. Let // be a constant matrix with 

tru = 0, [J, n] = 0, 

and associate to q in P a sequence of matrix- valued functions Fj. defined recursively 
by 

Fo{x) = n; 

dFu 

[J,Ffc+i] = — ^ + [Q,Ffc]; lim Fk+i{x) = Q. 
ax x^-oo 

The Fk depend nonlinearly on q for /c > 1 (/c > 2, if n = 2). Various formal and 



Theorem 2.1. [Sa, BC2, BC3]. Each Fk{q) is a polynomial in q and its deriva- 
tives of order less than k. The hierarchy of flows defined by 

(2.8) q = [J,Fk+i{q)] 

are Hamiltonian with respect to Vtp and the Hamiltonians are in involution with 
respect to the Poisson bracket { , }p. 

We shall discuss the Hamiltonians for these flows later. 

It is also well-known that the scattering transform linearizes the flows (2.8). We 
discuss this in the next section. 

3. The scattering transform; symplectic 
structure on scattering data 

We summarize here the basic results of scattering theory for the operator (2.1); 
cf [BCl]. In this section we assume 

(3.1) J = diag(iAi, .. .,iAn), Aj G R, Ai > A2 > ■ ■ ■ > A^. 

For a given q in P, we seek a matrix-valued solution of the spectral problem 

d 

(3.2) -—-ip{x,z)=zJtlj{x,z)-\-q{x)ip{x,z), zeC 
ox 

which is normalized by the asymptotic conditions 

(3.3) lim 'i/j{x,z)exp{—xzJ) = 1, linisup\\'ilj{x, z)exp{— xzJ)\\ < 00. 

If / [[^(x)!!^^ < 1, then there is a unique solution to (3.2), (3.3), and it has a limit 

(3.4) lim exp(-x^J)^(x, = ^(6, ^ G R. 

X — > + oo 

The transformation q \-^ s = s{-; q) is one of two versions of the scattering transform, 
and s is called the scattering matrix. 

Still assuming J \\q{x)\\dx < 1, the solutions i(j{x, z) for non-real z are holomor- 
phic and have limits on R which are related by 

(3.5) i;{x,^ + iO)=^|J{x,C-iO)v{0. ^ G R. 

To describe the target spaces for the maps q \-^ s and q \-^ v, we deflne the spaces 

(3.6) SL^ = {ae SL{n, C) : a^k = if ± {j - k) > 0}; 
SL^ = {ae SL^ : diag(a) = 1}; 

SL^ = {SL+ ■ SL-) n {SL- ■ SL+). 

This means that SL^ consists precisely of those s in SL{n) = SL{n, C) which have 
two (unique) triangular factorizations 

(3.7) s = s+vj^^ = S-vZ^, s± G SL , v± G SL^. 



(3.8) Definition. SD is the set of matrix-valued functions s : R ^ SL{n) with 
the properties 

(3.8a) s is smooth and hounded; each derivative has an asymptotic expansion in 

powers of ^~^ as \^\ -^ oo; 
(3.8b) s takes values in SL^, so it factors as s(^) = s±(Oi'±(0~"'^/ 
(3.8c) the diagonal-matrix-valued functions 5±(^) = diags±(^) are the boundary 

values of a diagonal-matrix-valued function which is hounded, holomorphic, 

and invertible in C\R. 

(3.9) Definition. SD' is the set of pairs of matrix-valued functions {vj^,V-),v± : 
R -^ SL^ , with the properties 

(3.9a) each entry of v± — 1 belongs to 5(R); 

(3.9b) the upper principal minors of the matrix-valued function v — vZ v+ are 
non-zero and have winding number zero. 

Condition (3.9b) excludes discrete scattering data (bound states); throughout 
this paper we consider only potentials with purely continuous scattering data. 
We equip SD and SD' with the Schwartz topologies. 

Theorem 3.1. (^[Sh], [B-Y], [BCl]j. The map q \-^ s is a diffeomorphism from a 
neighborhood of in P onto a neighborhood of 1 in SD, and it extends to map an 
open set in P bijectively to a dense open set in SD. 

The matrix function v in (3.5) is related to the matrix function s by the factor- 
izations (3. 7); in fact 

(3.10) V — f 1"*^^+ — sZ^Sj^, 

and this equation uniquely determines v± from v. The map q ^^ {vs^.V-) is a 
diffeomorphism from a neighborhood of the origin in P onto a neighborhood of 
(1, 1) in SD' , and it extends to map an open set in P bijectively onto a dense open 
set in SD' . 

(3.11) Remark. It is nearly implicit that SD and SD' are diffeomorphic. The fac- 
torizations (3.7) and (3.10) determine v from s. Conversely, write s± = 6±t± with 
5± diagonal and t± in SL^ . The factorization (3.10) gives t'Z^t'+ = tI^(5Z^5-|-)t+, 
showing that t± and 5Z 5-|- are determined algebraically from v. The holomorphy 
properties (3.9) show that the factors 5_, 5+ can be obtained from 5Z^5+ by solving 
a Riemann-Hilbert factorization problem, and thus s± and s itself can be obtained 
from V or from the pair {vj^, V-). 



Proposition 3.2. ^[G], [BCl];. The pull-back of the 2-form Op of (2.3) under 
the inverse of the scattering transform is 

(3.12) ns = ^ [ tr[v-{Sv)v^^ A s'^Ss]. 

It will be convenient to have a somewliat different formulation. 
Proposition 3.3. The 2-form O5 can be written 

(3.13) Vts = — I tr[v^'^5v+ A s+^5s+ - vZ^5v- A sl^5s_]. 



Proof. Since v = v_ V-\- and s = s±Vj. , it follows that 

v-{6v)v^'^ = {Sv+)v^'^ — {6v-)vZ^; 

s~^Ss = s~^[Ss±v^ — s±v^ Sv±v^ ] = v±s^ {5s±)v^ — {Sv±)v^ . 

Then 

tr[{5v±)v^^ A s~^5s] = tr[v^^5v± A s^^5s±], 

since {Sv±)v^ is strictly upper or lower triangular. 

Next we consider the image under the scattering transformation of the Poisson 
bracket { , }p of (2.7); equivalently, this is the Poisson bracket associated with the 
2-form O5 on scattering data. As usual, we may consider the entries of the scat- 
tering matrix s(^) = s(^; q) to be functionals on P and compute the corresponding 
bracket 

{sjk{0,sem{v)}si^) = {sjk{0,sem{v)}p{Q), ^,^e R-- 

There are two problems here. First, the gradients 5sjk{i)/5q do not decay, so the 
formula (2.7) does not have an absolutely convergent integrand and it is necessary 
to use a regularization such as 



lim rtrf[J,^(0]^(r7)Vx. 



Second, even this limit exists only in the sense of distributions in the two variables 
^, Tj. Thus the precise meaning of the calculation is this: for any pair of test 
functions u^w in C°°(R) one considers the pair of functionals 

Fiq) = / s,k{i)u{i)di, Giq) = / s,^{i)w{i)di, s = sUq). 



Then formally one has 

(3.14) 



SFJG. 



tr{[J,^]^)dx 
-N Sq'6q' 



{sjkiO^ siniir])}su{i)w{r])didr] 

as the defining equation for the distribution {sjk{i),S(,m{Ti)} ^ V^R. x R). The 
following calculation is standard; see [Ma] for the 3x3 case and [Sk], [KD] for 
i?-niatrix formulations. 

Proposition 3.4. The distribution defined by (3.14) '^^ given explicitly by 

(3.15) 

{sjkiO^^lTniv)} = T^iSjm{0sik{v)[s9n{i - j) - sgn{m - k)]5{^ - 7]) 

+ Sjk{0^im{v)[^je - ^km]p-v.——- 

where we take sgn{0) — and p.v. denotes the principal value. 
Proof. The variation of s with respect to q is 

(3.16) s{0 = f s{Oi^{x,0~'q{xmx,Odx; 

Jn 

[BC3, (2.45)]. Here the t/^ are the eigenfunctions (3.2), normalized at x = — oo. We 
write ifj^Xj^) = 'i/j{xj^)s{^)~^, which is normalized at a; = +oo. With F as above, 
an easy calculation using (2.2), (2.5), and (3.16) shows that 



SF f 

— (x)= ij{x,i)ekj^{x,i)~^u{i)di. 



A similar formula holds for G, so (3.14) becomes 

tr{[J.-^]-f)dx = 



lim /// tr{[J,tlj{x,^)ekj^{x,^)~^]ij{x,r])eme^{x,'r])~^)u{^)w{'r])dxd^d'r] 

N 



N^oo 




We use the identity 

1 d 



lil^ix^ri) ^'ib(x,f)]^'ib(x,r]) ^Jthix.Cl 



and the properties of the trace to conclude from (3.17) that 
(3.18) 

- ekjg{i, ?7, -N)emig{i, V^ -N)] 

in the sense of distributions, where g{^, rj, x) = '^{x, $,)~^'i/j{x, rf). Now 

g{i,r],x)^ey.Y>[x{r]-i)J]s{r]) as a; ^ +00; 
~ s(^) exp[a:(?7 — ^) J] as a; ^ — 00. 

Thus the right side of (3.18) is 

(3.19) lini ^[,^.^(^),,,(^)e^A^(^-'')(^^-^^) - .,>„(0«^fc(r7)e^^(«-^)(^--^'=)] 

There is no singularity in (3.19) since the term in brackets vanishes at ^ = ry; 
therefore we may replace the expression in (3.19) by the principal value integral, 
i.e. letting the distribution act as the limit as e | of the integral over the region 
1^ — ?7| > e. This allows us to decouple the two terms in (3.19) and use the identity 

lim p.i;.-^e^"^(^-^) = Tvisgn(a)d(C - v) if a e R\0. 

to deduce (3.15) from (3.19). 

As is well-known, the flows of Theorem 2.1 become linear on the scattering side. 

Proposition 3.5. The potential q{-, t) evolves according to (2.8) if and only if the 
scattering data evolve according to 

(3.20) _,(^,t) = e'=[/i,.(e,t)], g-^v±iU) = C''[t^,v±{C,t)]. 

For a proof, see for example [BC2], [BC3]. These flows are Hamiltonian (with 
the same Hamiltonian functions as in the original variables) with respect to the 
symplectic form O5 on scattering data since the structure has simply been trans- 
ferred from P to SD or SD' . We wish to emphasize that on the scattering side 
the flows (2.8) are not only linearized but decoupled for different values of ^, ry; 
equivalently, the Hamiltonian vector fields act in a pointwise fashion on the entries 
of s or of (z;+, V-). This allows us to reduce the question of complete integrability 
of the family of fiows (2.8)/(3.20) to a finite-dimensional problem. 

This integrability question is related to certain problems and questions concern- 
ing the symplectic and Poisson structures on scattering data. Observe that the 
2-form O5 lifts to the loop space 
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More precisely O5 is the puUback to SD of the 2-form defined by (3.12) or (3.13) 
on L{SL^). Moreover, these are pointwise formulas, in the obvious sense: they 
express the form as the direct integral of forms computed pointwise from entries of 
a in L{SL^). However, the form is not symplectic on L{SL^). (This can be seen 
from the fact that the dimension of SL is n^ — 1 but the rank of the pointwise form 
in (3.13), as we show in the next section, is n^ — n.) 

The alternative space of scattering data, SD' , is itself a loop space, whose point- 
wise dimension is n^ — n, and the form O5 is symplectic on SD'; however the 
computation of s from {v+,V-) involves the solution of a Riemann-Hilbert prob- 
lem, so that (3.12) does not express O5 as a direct integral of pointwise 2-forms on 
the fiber of the loop space; in fact it has no such pointwise expression when n is 
larger than 2. 

The Poisson bracket { , }s lifts to the whole loop group L{SL{n)) using the same 
formula (3.15), but the nonlocal term involving p.v.{C '~ ^)~^ shows that { , }s is 
also not a direct integral of a pointwise defined bracket, unlike the bracket { , }p 
on the space P of potentials. This corresponds to the fact that the submanifold 
SD C L(SL{n)) is determined in part by the nonlocal constraint (3.8c). The 
nonlocal term does not vanish when { , }§ is considered as a bracket on SD' if 
n > 2, which proves the earlier assertion that fls has no pointwise expression on 
SD'. 

The following summarizes these observations. 

Proposition 3.6. O5 is a dosed, pointwise 2-form but not a symplectic 2-form on 
the loop space L{SL^). O5 is a symplectic 2-form hut not a pointwise 2-form on 
SD'. 

{ , }s is not a pointwise Poisson bracket on L{SL{n)), nor on SD, nor on SD' 
(when n > 2). 

Our main positive results involve choosing new coordinates to overcome the 
limitations in Proposition 3.6. To state them we must extend the notion of a 
distribution-valued Poisson bracket beyond the coordinate functions themselves. 
Suppose /, (7 are in C°^{SL{n)) and u,w are test functions. We define functionals 
on SD, or on P, by 

/?(«) = /(«(0), F{s) = lf{s{OHOdC. 

9d^) = 9{s{0), G{s) = l9{s{0)w{0dt 
Thus the coordinate functions ajk{s) = Sju, a£m{s) = S£m give rise to the func- 
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distribution {/^, (7^} G 'D'(R x R) by the formal calculation 

(F,G)p = JJ{k.gr,}su{^)wirj)dCdrj. 

This leads also to the expression 

Ubjl^ Ubfi^jyi 

Theorem 3.7. There are functions p,^,qi,,l < u < \{n'^ — n) which are defined 
and holomorphic on a dense open algebraic subset of SL{n) and which have the 
following properties. 

(3.21) The map s 1— ^ (pi o s, p2 o s, • • • , Qi o s, ^2 o s, • • • ) is an injection from a dense 
open set in SD into C°°(R; C'""). 

(3.22) Os = E.5(p,os)A5(g,os) 

(3.23) {p^,^,Piy,y}s = {qiy,^, qu,ri}s = 0; {Pf,,^, Q^,^}s = 6f,J{C - 77) . 

(3.24) For any f in C*{SL{n)), the distributions {Pv,(^, fn} have support at ^ = r], 
all V . 

In other words, the functions p,y o s,q,y o s are global Darboux coordinates on 
the manifold SD of scattering data. The additional fact (3.24) implies that Hamil- 
tonians which are functions of the Pi, give pointwise vector fields on SD. These 
same functions Pi, provide a strong positive answer to our question about complete 
integrability in the SL{n) case, as follows. 

Theorem 3.8. The functions p^, of Theorem 3.2 may be chosen so that for each 
traceless diagonal matrix fx, the Hamiltonian for the flow (2.8), (3.20) is a linear 
combination of the functional 

(3.25) [ epAs{0)d^- 

Thus the family of flows (3.20), which is determined pointwise by the (n — 1)- 
parameter family of traceless diagonal matrices, is imbedded in the |(n^ — n)- 
parameter family of flows generated by the p^^s. In fact the functions p^jQu of 
Theorem 3.2 provide action-angle variables for the flows (3.20). 

Theorems 3.7 and 3.8 are proved in section 6. 

In this section we have made two restrictive assumptions about J : that the 
eigenvalues are distinct and that they lie on a line through the origin in C. The 
same results hold without the second assumption, as we show in section 8. 

Another situation arises with reductions, i.e. the imposition of restrictions on the 



12 

q + q* =0, still assuming (3.1). For such q, the scattering data satisfies the 
corresponding constraints 

s(0 belongs to SU{n); v+{0*V-{0 = 1. 

The pullbacks of the 2-forms are still symplectic. When n = 3 the simplest associ- 
ated nonlinear evolution equation is the 3-wave interaction. (Note that one needs 
fx + fx* = in (2.8), (3.20) to preserve the constraints.) Here our manifolds are real, 
and we need the functions p,^, q,y of Theorem 3.2 to be real in order to preserve the 
structure. 

Theorem 3.9. The functions pj^^q^ of Theorems 3.7, 3.8 can he chosen to he real 
on SU{3) in the case n = 3. In particular, the three-wave interaction is a completely 
integrahle Hamiltonian evolution in the strong sense. 

This result is proved in section 7. 

4. A 2-FORM IN GL{n) 

Let GL{n) denote either GL{n^ R) or GL{n, C); in the latter case the functions 
and forms to be considered are complex- valued. The key steps in deriving Theorems 
3.2, 3.3, and 3.4 involve an analysis of a 2-form in GL{n) : the integrand in (3.13). 
As in section 3 we introduce the matrix spaces 

GL^ = {ae GL{n) : ajk = if ± {j - k) > 0}; 

GLJ = {ae GL^ : diag(a) = 1}; 

GL^ = {GL+ ■ GL-) n {GL- ■ GL+). 

Again, GL^ consists of those elements of GL{n) with factorizations 

(4.1) a = a+v^'^ = a-vZ^j a± E GL , v± G GLj. 

We consider a±, v± here as functions of a in GL{n); then da± and dv± are matrices 
of 1-forms on GL{n) and we may define a 2-form by 

(4.2) O = tr[v^ dv+ A a^ da+ — vZ dv- A aZ da-]. 

At the identity, = 2^ dajk A dakj, so O has rank > n^ — n on an open set. 

j<k 

Theorem 4.1. On a dense open algehraic suhset of GL^ the 2-form O has a rep- 
resentation 

^ 1 

(4.3) n = ^dp^hdq^, N=-{n^-n), 

where p^, and q^, are analytic (holomorphic) and the 1-forms dp^^dq^, are indepen- 
dent. In particular, O is closed and generically has rank v? — n. 

The proof of this theorem is given after Theorem 4.4. The strategy is to obtain 
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Lemma 4.2. For n — 2, 



(4.4) 



fl — d log 



011022 

A 



A dlog 



021 

ai2 



A — det a. 



Proof. If o is in GL^{2) then 



^+ 



1 

-021/022 1 

1 -012/011 
1 



A direct calculation gives 



n = d 



-d 



'a2i' 
.022. 
'ai2 

_Oii_ 


A 
A 



_ A/022 ai2 
0+ — 

022 

Oil 

021 A/oii 



oi2(io22 a22dai2 



A 
0210011 

A 



A 

0110021 

A 



and some further manipulation leads to (4.4). 

(4-5) Remark. Given I < j < k < n, consider the subset of GL* consisting of those 
o whose only non-vanishing off-diagonal entries occur in the (j, k) and (/c, j) places. 
The puUback of O to this subset has, by an analogous computation, the form 



(4.6) 



O = dlog 



O'jjCtkk 



A dlog 



lO'jk^ 



A 



UjjUkk ^jk^kj 



We now proceed to reduce the general case to a sum of cases as described in the 
preceding remark. 

Suppose that tti, 7r2, . . . , tt^t, A^ = ^(n^ — n), are permutation matrices with the 
two properties (written using the standard matrix units) 

(4.7a) TVi, is the matrix of the transposition (k,k + 1) , k = k^; 

(4.7b) the product 7ri7r2 . . . tvn is the antidiagonal matrix r = ^ Cj^n+i-j- 

There are various such decompositions of the antidiagonal matrix r. One such 
decomposition corresponds to permuting (1,2, . . .,n) by moving 1 to the extreme 
right in n — 1 steps, then moving 2 to the position left of 1 in n — 2 steps, and so 
on. 

(4.9) Definition. Given n^, satisfying (4 ■'7), set 



ro = 1, ri, = 7ii7i2- --T^u, l<i^<N; 

TT „ /^r+„ — 1. r . /^r— „ — 1. 7-1 
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where B^, = {b E GL{n); bjk = for j ^ k unless {j, k} = {k^,, k^, + 1)}. Thus the 
matrices in D^, are block diagonal in the sense that after conjugation by r~^ the 
nonzero entries lie in a single 2x2 block along the diagonal. 

Note that U^, L^, D,^, U^+D,^, and Li^+D^, are subalgebras of the matrix algebra 
M^. Let 

(4.9) Pu-.Mr,^ D, 

be the projection which commutes with left and right multiplication by diagonal 
matrices. Then 

(4.10) {D, + U^) n {D, + L,) = D^; 

(4.11) Pj^ is an algebra homomorphism on D^ + U^, and on D^ + L^^. 

Note also that because r^-i and ri, differ by a single transposition, we have the 
identities 

(4.12) D, + L,_i =D, + L,, D, + t/,_i = D^ + U,. 

(4.13) Definition. Suppose the permutation matrices Wi, satisfy (4-7), and suppose 
r~^ari, is in GL^, < u < N. Then we may factor the r~^ari, as in (4-1) io obtain 
matrices u^,£i, such that 

(4.13a) aUi, = i^y, u^, G U,y, ijy E L,^. 

Now define v~^ ,v~ , a+, a~ and 2-forms O,^ as follows: 

(4.13b) 

v~ = P^u^-u vl = PuU^, a~ = PJu-u a'^ = PJu] 
(4.13c) 

^v = tf[{'^t)~^dv^ A {a'^)~^da'^ — {v~)~^dv~ A {a~)~^ A da~]. 

Lemma 4.3. With a as in (4- 13), 

(4.14.) V- = Wo, a- = £oj v+ = un, «+ = ^n 



(4.15) uZ^_^u, = K\i. = {v^r^vt = {a-)-^a 



1^+ 



In particular, a+(f j") ^ and a^ {v^ ) ^ have a common value a^, in D^, and Oj^ is 
the pullback at a^, of Vt under the map b ^^ r~^bri, from D^ to B^ C GL{n). 

Proof. That (4.14) holds is clear from the definitions, together with the assumption 
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(4.15), since i,y-iu~_i = a = lj,u~^ . Because of this first identity and (4.12), the 
common value belongs to {Dj, + L^j) fl {D^, + Uij) = D^. Therefore we can project 
and use the property (4.11) to obtain the remaining identities in (4.15). The final 
statement is immediate from conjugation of Q,y by r,^. 



Note that a,^t'jf = aj^, and that Q^, can be expressed directly in terms of the 
entries of the tti, via (4.6). By virtue of the following decomposition theorem, the 
computation of O is reduced to a sum of 2 x 2 problems, as in Remark 4.5. An 
algorithm for computing a^, in terms of a will be given below. 

Theorem 4.4. Under the assumption (4-V' ^ ^-^ ^^e sum 

(4.16) = 0l + 02 + --- + 0Ar. 

Proof. We begin by reversing the reasoning in the proof of Proposition 3.3 to write 
O in the alternative form 

fl = tr[v-{dv)v^ Aa~^da]. 

From (4.14), (4.15) it follows that 

V = VZ V+ = Uq Un = 1'l'l'2 • • • vn, 

where f j, is the common value of the matrices in (4.15). Then 

v-{dv)v+ = v-{^{vi . . .v^-i)dv„{v„+i ...vn)}v:^'^ = ^u„-i{dv„)u~'^. 

Note the identity 

u,,-idv^u~^ = u^{v^)~^dv^u~^ 

-u^-i{v~)~'^dv~u~l^. 

From a = £i,u~^ we find 

tr{uiy{v'^)~^dv'^u~^ A a~^da} 

= tr{{v+)-^dv+ A [{du-^)u^ + i-^di^]} 

The first term vanishes since v^ and u^, both belong to Uiy. For the second term 
we have 

tr{{v+)-^dv+AP^{i-^di,)} 
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since P^, is multiplicative on L^,. The second term in tr{u,y-idv^Uj^^ A a ^da) is 
treated in the same way, using a = ijy-iu~_i, and (4.13) follows immediately. 

Proof of Theorem 4-P Choose permutation matrices which lead to a decomposition 
(4.16) of O. According to Lemma 4.2 and Remark 4.5, each Q^, has the form 
dpi, A dq,y, so we obtain the desired representation (4.3) on the dense set where the 
Pv,qu are defined. It follows immediately that Vt is closed and that it has rank at 
most ri^ — n everywhere. The rank is n^ — n at the identity, and therefore is n^ — n 
on a dense algebraic open set, so the dpj,^ dq^, are generically independent. 

To compute Vt^, and hence p^, q^, we need to find a,y. To obtain a^, from a, note 
that 



(4.17) 
where 
(4.18) 
and 



au 



# 



?# 
'i^ 



u 



# 



Uu-l{v^) \ if =i^-i{v^) \ 



PJ* = PJ^-i{v^) ^=a^{v^) ^=a^, P^uf = l. 



The matrix r~_j^a,yr,y_i is block diagonal; its nontrivial part is the 2x2 block with 
entries from rows and columns ki,, /cj^+i . It follows from (4.17), (4.18), and a simple 
computation that this block is D — CA~^B, where A is {kv — 1) x {ki, — 1), D is 
2x2, and 

r -1 . f A B 



r^_^ari,_i\j^k<K+i - \ (j jj 

The following notational convention will be useful. 

(4.19) Definition. If J and K are two subsets of {1,2, ... ,n} having the same 
cardinality and a is in M^, then m{J] K) = m{J; K; a) denotes the determinant of 
the corresponding suhmatrix 

m{J;K;a) = det{ajk)jeJ,keK- 
We seim(0;0;a) = 1. 



Direct calculation leads to 

1 



D-CA-^B 



■m{J,k^]J,kj,) m{J,k^] J,k^ + 1) 

m{J,k„ + l;J,k^) m{J,k^ + l;j,k^ + 1) 



m{J; J) 
where J = (1, 2, . . . , /c,^ — 1) and m{K; L) = m{K; L; r~_]^ar;^_i). Thus 

O,, = dpi, A dqi, 

m{J, ku] J, ku)m{J, k^, + l;J,k,y + l 



(4.20) 



Pi^ = log 
q„ = log 



_m{J; J)m{J,kt,, k^ + 1; J, k^^k^, + 1)_ 
m{J,ki, + 1; J, k^y 
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(4-21) Example. For n = 3, we consider the decomposition of r corresponding to 
(123) -^ (312) -^ (231) -^ (321), i.e. 



"0 1 0" 




'1 0" 




"0 1 0" 




'0 1" 


1 




1 




1 


= 


1 


1 




1 




1 




1 



Let A be the cofactor matrix A = A(a ^)*, where A = deta. Then the corre- 
sponding decomposition of O is 

O = d\og{— ) A dlog{ ) 

, ,, /^11^33n . ,, /^13^ 

+ dlog(— )Arflog(- 



+ c^log( 



Aa22 



■^31- 

Arflog( — ). 

^23 



We conclude with some symmetry properties of O which wiU be important later. 
Proposition 4.7. The 2-form O is odd under the automorphisms of GL(n) 

$i(a) = {a~^)\ $2(0) = rar, 

(where r is again the antidiagonal matrix '^e.j^n+i-j), i-e. 

$*0 = -0, J = 1,2. 

Proof. ^j{GL^) = GL"^. Therefore if a^ = $j(a), the ± factors in (4.1) are $j(a=p), 
$j(i'=p). It is immediate from this that $2^ = — O. The result for $1 makes use 
also of the identity ^l{b~^db) = —{b~^dbY, together with the identity tr{a A /3) = 
tr{a^ A /3*) for matrix-valued 1-forms. 

5. A SYMPLECTIC FOLIATION AND POISSON BRACKET ON GL{n); FLOWS 

We introduce now a foliation of GL^ which is naturally associated to the factor- 
izations (4.1). As in remark (3.11) we define diagonal matrices d± = diag(a-i-) and 
S = 5Z^5_|_, and set 



(5.1) 
so that 
(5.2) 



a± = S±b±, b± G GLq 



vZ^v+ = aZ^a+ = bz\6Z^6+)b+ = bZ^6b+. 



In the notation of (4.20) the principal minors of a G GL{n) are m{J; J]a)] we 
abbreviate this to m{J;a). If a is understood, we may write m{J). In particular 
the upper and lower principal minors are 



/'K 0\ 



J+ I+/„\ 



,lt 



A\. 



J A 



.\. 
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It follows from (4.1) that 5_|_(5_) has the same lower (upper) minors as a, so 

(5.4) (S+),, = -^; (5.),, = -^, l<j<n. 

Therefore 5+ and 5- are determined uniquely by 5 = 6Z^6+, together with the 
quotients 

(5.5) (pj{a) = d+{a)/dj^^{a), I < j < n. 

Note that (/'n(tt) = deta, and that the Lpj determine the product S+S-. 

Theorem 5.1. The foliation of GL^ by the functions (^j is symplectic for Vt . Each 
leaf {a : (pj{a) = Cj, 1 < j < n} is parametrized by F* = {(f^+, V-) E GLq xGLq : 
vZ^v+ E GL}; the pullback of fl to a leaf is generically of rank n^ — n; the pullback 
from a leaf to F* is independent of the choice of leaf and is given by 

(5.6) tr[v+^dv+ A {b+^db+ + b^\d-^d6)b} - vZ^dv- A {bZ^db- - bZ\6-Ud)b-)}] 
where b± and d are determined from (v^jV-) by the factorization 

(5.7) vZ^v+ = bZ^Sb+, b± E GL^, S = diag(5). 

Proof. Starting with a in GL*, we define 6±, 5±, as above and again set 5 = 5Z^5-\-. 
Let rj = d-d+. Then 

6^ d6± = ri~^dri ± 6~^d6] a^ da± = b^ db± + b^ {6^ d6±)b±. 

Consequently O is given by the sum of (5.6) and 

(5.8) tr[v^^dv+ A b'^^{r]~'^dr])b+ - vZ^du- A bZ^{v~^drj)b-]. 

The pullback of dr] to the leaves of the foliation determined by the functions (fj 
vanishes, since these functions determine ry, so the pullback of O to the leaf is given 
by (5.6). Now b± and 6 are determined from (f_|_, t;_) in V^ by (5.2), so the pullback 
of O to y* is leaf-independent and given also by (5.6). These puUbacks have rank 
< dim(KK) = n^ — n everywhere and rank n? —n dX the unique diagonal element in 
a given leaf, so they have rank n? — n generically. 

The symplectic foliation gives rise to a Poisson structure on GL^. In fact there 
is a Poisson bracket ( , )l on each leaf L corresponding to the pullback of O to L, 
and this may be extended to a (degenerate) Poisson bracket for functions on GL^, 
characterized by 
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Equivalently, if fi == T,dpi, A dq^, as in Theorem 4.1, then 

(5.10) (Pf^.qu) = S^:,, {Pii:Pu) = Q = {qi,,qu), 1 < /«, ^^ < (n^ - n); 

(/,V?,) = 0, l<j<n. 

Functions such as the (pj which Poisson commute with aU functions are sometimes 
called Casimirs. 

This Poisson bracket was computed for the standard coordinates of GL{n) by 
Lu [Lu] in the cases n = 2, n = 3; Lu also conjectured the general form below 
and pointed out the connection with the classical limit of a quantum version due 
to Drinfeld [Dr], as noted in section 1. 

Proposition 5.2. The Poisson bracket (5.9) is odd under the automorphisms 
^1,^2 of Proposition 4-7, i.e. 

(/o$„ (70 $,) = -(/, (7)0$,, J = 1,2. 

Proof. The Casimirs (fu satisfy (fk o $j = fn+i-k, so $j maps leaves to leaves. 
Proposition 4.7 implies that the puUback under $j to a leaf L of O on ^j{L) is — O. 
Therefore the pushforward of the Poisson bracket ( , )l is — ( , )^^{L)^ and (5.9) 
implies the desired result. 

Theorem 5.3. The Poisson bracket on GL^, induced by the 2-form O and the 
foliation by functions {(fj) extends to the full matrix space Mn and is given by the 
following bracket relations between the coordinate functions ajk'. 

(5.11) {ajk, aim) = j[sgn{(.-j) - sgn{m - k)]ajma£k, 

where again sgn{0) = 0. 

Proof. We show first that the calculation can be reduced to the cases n < 4. For 
a fixed k, 1 < k < n, consider the map Mn -^ M„_i obtained by omitting the 
K-th row and column. We claim that the pushforward to M^-i of the Poisson 
structure on M^ coincides with the structure on M^-i itself, i.e. the Poisson 
bracket of coordinate functions ajk, aim on M^, i, j, /c, £ 7^ k, is the same as that 
obtained by considering them as functions on M^-i- To verify this claim we take the 
decomposition (4.7) of the antidiagonal matrix r G M„ obtained from the following 
three sets of permutations: the first set takes (l,2,...,n) to (1,...,k— 1,k + 
1, . . . , n) in n-steps; the second set takes us to (n, n — 1, . . . , 1, k) in \{n — l)(?i — 2) 
steps; the third set takes us to (n, n — 1, . . . , 1) in [k — 1) steps. The corresponding 
additive decomposition of O then takes the form 
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with the obvious notational convention. According to the prescription after Lemma 
4.5, O" = Tidpi, A dqi,, where the pi,, q^, are functions of the matrix h in Mn-i which 
corresponds to a in Mn by the map above; moreover O" has exactly the same form 
as 0^_i. To complete the verification of our claim it is, therefore, sufficient to 
show that the foliation functions for h Poisson commute with all entries of h when 
these are considered as functions of a. We know that the p^^ which correspond to 
O' and to Q." in the decomposition (5.12) commute with all Pv,Qv from O" and 
with each other, so it is enough to show that the foliation functions of h in M^_i 
are computable from the pi, in O' and O'", together with the foliation functions on 
Mn- The Pj, corresponding to O' are 

log(fi'fc+i/fi'fc), K<k<n; fiffc =m(l, ...,K, ...,/c)/m(l,...,/c) 

where again m{. . .) denotes the principal minor of a based on the indicated rows 
and columns. Similarly, the pi, associated to O'" are 

log{hk+i/hk), 1 < k < k; hk = m{k, . . . ,n)/m{k, . . . ,k, . . . ,n). 

Modulo the foliation functions on M^, the gk and hk can be determined from the 
p^. Again modulo the foliation functions on Mn, the gkS and h^s are equivalent 
to the set of functions 

771(1, ..., K, ... , k)/m{k + 1, . . . ,n), k < k < n; 
m(l, . . . , A; — l)/m(/c, ...,«:,... ,n), 1 < k < k. 

These are precisely the foliation functions of b as an element of M^-i- This com- 
pletes the proof that the two Poisson structures coincide on b. 

Suppose now that aij^aki are any two coordinate functions on M^. Repeated 
use of the preceding argument shows that their Poisson bracket can be computed 
by taking them to be functions on Mp, with p the cardinality of {i,j, k,£}. Thus 
the computation is reduced to the cases Mi (trivial), M2,M3,M4. The complete 
computation is tedious, so we merely indicate a few representative cases. Recall 
that a Poisson bracket is a derivation for each of its arguments. 

For 71 = 2 the foliation functions are aii/a22, A = deta. From this fact and 
(4.4) we deduce 

(aii,a22) = (aii,aii(a22/aii)) = {0,11,(111)0-22/0-11 = 0; 
(an, ai2a2i) = (an, ana22 - A) = 0. 

Also, p = log(ana22/A), q = log(a2i/ai2) and {p, q) = 1, so 

(^11^227^21/^12) =^11^22^21/^12; 

(an,a2i/ai2) = (a^^, a2i/ai2)/2an = (011022^11/^22, a2i/ai2)/2an 
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Therefore 

1 1 

(aii,a2i) = 7" (aii,a2i = 012(^115 tt2i/fli2)/2 = -011021 

2a2i 4 

Because of the symmetries in Proposition 4.7, (022, «2i) = —(011,^21) and so on. 
Also 

(021,012) = (021, ai2a2i)/a2i = (a2i, — aiia22)/a2i 

— (^11^22^11/^22, «2i)fl2i/fliifli2 = 2(aii, 021)022/^21 = 011022/2. 

For n = 3, reduction to n = 2 gives all brackets such as (022, ^12)7 (^237 ^33)- Let 
A = (deto)(o~^)* be the cofactor matrix. The foliation functions are 

deto. On /All, 033/^33 • 
The decomposition (4.21) implies, therefore, that 

= (oi2aiiOi3/AiiA33, A13/A31) = 2o22aii033(o22, ^13/^3^/^11^33 

(^22, Al3^3l) = (0^22, ^11^33 - 0^22 det o) = ^11^33(022, aii033)/oii033 = 0. 

Therefore (022,^12) = = (022,^31), and the 2x2 results allow one to calculate 
(o22,tt3i) and (022,^13) from these identities. Similar computations yield all the 
3x3 brackets, though for some it is convenient to replace the decomposition in 
(4.21) with the decomposition obtained from factoring r by means of (123) -^ 
(132) -^ (312) -^ (321). 

The case n = 4 is similar, and again brackets like (012, 024) are known from the 
3x3 computation. This completes our sketch of the proof of Theorem 5.3. 

We consider now the (n — 1) -parameter family of flows in M^: 

(5.13) o(t) = exp(t(u)o(0) exp(— t|u), fx diagonal, trn = 0. 

This conjugation preserves the principal minors of o, so the flow preserves GL{n), 
GL^, and the leaves of the foliation. The factorization (4.1) is also preserved, so 
the flow preserves the 2-form O. Therefore these flows are Hamiltonian. 

Theorem 5.4. The Hamiltonian function for the flow (5.13) is tr[ii\ogS\, d = 
6Z S-\-, and it is a linear combination of the functions p^, of (4-20). 

Proof. We use the additive decomposition of Theorem 4.4, with the factorization 
of r described after (4.7). If n^ is associated with the interchange of j and /c, then 
(4.20) implies that under the flow (5.13), Qu = Hk — fJ'j- Therefore the Hamiltonian 
for (5.13) is a linear combination of the pi,. The pi, themselves are logarithms of 
quotients of principal minors of o: 
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The term rrijk occurs in the numerator of exp(p;^) when yr,^ is associated to the 
interchange of j — 1 with k or of j with k+1; the term rrijk occurs in the denominator 
when n^ is associated to the interchange of j — 1 with k + 1 or of j with k. Thus 
the total weight attached to log mjfc in the Hamiltonian for the flow (5.13) is 

with the convention that if j — 1 = or /c + 1 = n, the corresponding term in 
parentheses is omitted. Thus logrujk has weight zero unless j = 1 or k = n, and 
the Hamiltonian for (5.13) is 

n— 1 n— 1 

"^{Hk+i - IJ-k) logmiA; + ^(/Uj+i - iJj) logrrijn 
k=i j=i 

n n 

= Xl^J'^°s[t?7(a)c?+_i(a)/c?7+i(a)rf+(a)] = ^^jlog%. 
i=i i=i 

(5.14) Remark. Theorem 5.2 and its proof show that the flows (5.13) are completely 
integrable in the classical sense: they are an (n — l)-parameter family of commuting 
Hamiltonian flows in each n? — n dimensional symplectic leaf of the foliation, and 
are part of the \{n^ — n)-parameter family of commuting flows generated by the 
Pj^. Note that the flows (5.13) are the only members of the larger family which are 
linear as flows on the full matrix algebra M„; in fact the generator of a linear flow 
which commutes with all the flows (5.13) must have each matrix unit Cjfc, j 7^ k, as 
an eigenvector and if such a flow leaves all the p^, invariant it can be shown to be 
included among the flows (5.13). 

6. Proofs of Theorems 3.7 and 3.8: 
Darboux coordinates for scattering data 

Up to a trivial normalization, the functions pi,, Qi, of Theorem 4.1 are the func- 
tions of Theorem 3.7 and Theorem 3.8. To see this, we return to the notation 
introduced before Theorem 3.7. Observe that if /, (7 belong to C{SL{n)) then in 
view of Proposition 3.4 the distribution- valued Poisson bracket can be decomposed 
as 

1 



(6.1) {/^, g^]s = [/, g]{sm5{i - 77) + (/, 9){s{i),s{r^))p.v.-^ _ 

where [ , ] and ( , ) have the following properties. 
(6.2) The map f,g ^-^ [f^d] is an alternating bilinear map from C°°{SL{n)) x 
C°°{SL{n)) to C^{SL{n)) which is a derivation in each variable: 

[fg,h]=f[g,h]+g[f,h] 
r -<• „i,i „r-f 1,1 I j,r-f „i 
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(6.3) [ajk, aim] = Tviajmaik[sgn{£ - j) - sgn{m - k)]. 

(6.4) The map f,g ^ (1^9) is a symmetric bilinear map from C^{SL{n)) x 
C^{SL{n)) to C^iSLin) x SL{n)) such that 

{fg, h) {s, s') = f{s) {g, h) {s, s') + g{s) (/, g){s, s') 
(/, gh){s, s') = g{s'){f, h){s, s') + h{s'){f, g){s, s'). 

(6.5) {ajk,aem){s,s') = ajk{s)aem{s')[Sje - Skm]- 

Here the ajk are the coordinate functions and s, s' are points of SL{n). 

The properties (6.2), (6.3) imply that the bracket [ , ] is precisely 47rz( , ), where 
( , ) is the Poisson bracket (5.11); this corresponds to the fact that O5 is the direct 
integral of AniQ.. Consequently we may replace the functions pj^ and q^^ of section 
4 by the renormalized versions 

(6.6) -p^, —-q^ 

to obtain 

(6-7) {Vp.,i, qv,r^}s = S^J{^ -V) + {Pt^, Q<y){s{0^ s{ri))p.v.-—— 

To complete the proof of Theorem 3.7 we must show 

(6.8) (p.,/) = 0, all / in C^{SL{n)), aU v- 

(6.9) (<?M'?'^) = 0' ^11 ^^■:^■ 

(6.10) Definition. Given subsets J, J', K, K' of {1,2,..., n}, set 

e{J, K- J', K') = card{J f] J') - card{K f] K'). 

Lemma 6.1. The bracket ( , ) between minors of s satisfies 

(6.11) {m{J; K), m{J'; K')){s, s') = e{J, K; J', K')m{J; K; s)m{J'; K'; s') 

Proof. The case when J, K, J', K' all have cardinality 1 is immediate from 6.5. 
The general case follows by expanding the determinants and using the derivation 
property 6.4. 

We can now prove (6.8) and (6.9), using (4.20). Each p^, is the logarithm of a 
term m{J)m{J') /m{K)m{K') , where each element j of (1, 2, ..., n) occurs with the 
same frequency in the pair of sets J, J' as in the pair K, K' . From the derivation 
property (6.4) we deduce that (6.11) is equivalent to 
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Each coordinate function ajk is itself a minor, and therefore (6.12) imphes 

m{J)m{J') 
log logajfc 

m[K)m[K') 

= e{J,J;j,k)+eiJ', J'; j, k) - e{K, K; j, k) -e{K',K'- j, k) 
= 0. 

This proves (6.8). To prove (6.9) we note that according to (4.20), each q^, has 
the form logm( J; K)/m{K] J), which we abbreviate slightly as log m{JK)/m{K J). 
Again we deduce from (6.12) that 

mjJK) mjJ'K') ^ 
log^— — , log 

m[KJ) m[K'J') 

= e{J, K- J', K') - e{J, K; K', J) - e{K, J; J', K') + e{K, J; K\ J') 
= 

because e{J,K;J'K') = -e{K, J; K', J'), and e{J, K; K', J') = -e{K, J; J', K'). 
This proves (6.9). For the injectivity property (3.21) we need the foliation functions 
ipj in addition to pj^, Qjy. By Theorem 5.4, entries of 5Z^5_|_ are linear combinations 
of the Piy. As in Remark 3.11, 5- and 5+ are Riemann-Hilbert factors of 5Z 5_|_. 
Finally, (5.4) and (5.5) determine the (fj from 6j^,d-. D 

Proof of Theorem 3.8. It is a well-known fact that the Hamiltonian for the flow 
(2.8) is the negative of the coefficient of z~^~^ in the asymptotic expansion of 

triJL\og6{z) as 2 ^ oo, where 

5{z) = lim 'i/'(3^7 -2) exp(— x^J), ^ G C\R. 

X — » + oo 

See [S], [BC3]. Now 5 is piecewise holomorphic with limits 5± on R, so the Hamil- 
tonian can be expressed in the form 



-7^ I etr^i\og[6z\i)5+{i)]di. 
2txi Jr 



According to Theorem 5.4, this integral is a linear combination of the integrals 

(3.25). 

7. Coordinates and flows on 5't/(3) 

Formula (4.2) defines a complex 2-form on the intersection of GL = GL{n, C) 
with the real submanifold GL^, fl U{n). The automorphism a 1-^ (a~^)* takes GL^ 
to GL^. Therefore, in the factorizations (4.1), we have 

(7.1) a± = (a-^)*, v± = {v-^y, aeU{n). 

A „ ;„ j-U„ c „f T3 „;j-; — A >-7 „.„ „„„ „ .1,, J„ j-U„j- r> T^ — r~<T r\ tH ^\ 
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Theorem 7.1. (a) iVt is a closed real 2-form on GL^ fl U{n). 

(h) The foliation of GL^ fl U{n) induced by the foliation of GL^ in section 5 has 
leaves with real dimension v? — n. The 2-form generically has rank v? —n on each 
leaf, so the foliation is again symplectic. 

(c) If n is a real diagonal matrix with tr{n) = 0, the flow 

(7.2) a{t) = e**^a(0)e-^*^ 

is Hamiltonian in GL^ fl U{n) with real Hamiltonian function tr{nlogd), where 
again 6 = 62^6+, 6± = diag(a±). 

Proof. Part (a) follows from the preceding remarks. For part (b), note that (7.1) 
implies 

(7.3) 6±{a) = [6^{a)-^]*, a G U{n). 

Because of (5.3) and (5.4), (7.3) implies that the foliation functions (5.5) take values 
in {l^l = 1}, so the induced foliation is defined by the n independent real functions 
arg ipj . The 2-form iQ has (real) rank n^ — n at the unique diagonal element in 
each leaf, so it has rank n? — n generically. Finally, the Poisson bracket associated 
to this symplectic foliation is — z( , ), where ( , ) is the Poisson bracket of section 
5 restricted to U{n). According to Theorem 5.3, therefore, the Hamiltonian for the 
flow (7.2) is tr(/ulog5), and according to (7.3) S = SZ 5+ = 5!j_5-|- is real. 

The Darboux coordinates Pi,, q^, constructed in section 4 are not real when spe- 
cialized to U{n) (and suitably normalized) except for n = 2. We show in this section 
that real Darboux coordinates can be chosen in SL{3, C) in such a way that: (a) 
the restrictions to SU{3) are real; (b) linear combinations of the p,y still include the 
Hamiltonians for the flows (7.1). As we shall see, the third Hamiltonian flow which 
commutes with the 2-parameter family of linear flows (7.2) is not linear on M3. 

We begin our discussion by recalling the Darboux coordinates for SL{3) in ex- 
ample (4.21). Corresponding to the new normalization zO, we take these to be 

pi = log(aiia22M33), p2 = log{AiiAss/a22), Ps = log(a22a33Mii); 
gi = ilog(a2i/ai2), i?2 = ilog{Ais/Asi), qs = ilog(a32/a23)- 

Here again the Ajk are the entries of the cofactor matrix (deta)(a~^)*. Therefore 
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It is convenient to make a preliminary linear canonical transformation to new Dar- 
boux coordinates 

(7.6) pi =pi +P2 =log(aiiAii) 

P2=Pl+P2+P3 = log(aiia22a33) 

P3=P2+P3= log(a33A33) 

qi= q2-q3 =i log(a23^13/a32^3l) 

q2 = qi-q2+q3 = «log(a2lA3ia32/ai2Ai3a23) 

qs = -qi + q2 = zlog(Ai3ai2M3ia2i). 

Then pi and ps are real on SU{3), but p2 is not. It follows either from direct cal- 
culation using the Poisson bracket (5.11) or from the first symmetry in Proposition 
5.2 that log(AiiA22^33) is also in involution with pi and p2, so we may take 

(7.7) Ij^ttjjAjj, j = 1,2,3, 

as the action variables for a new set of Darboux coordinates. The corresponding 
angle variables Gi,G2,©3 are then obtained by the classical Liouville method; cf 
[W]. We briefiy recall the method: in principle one solves the equations 

(7.8) pj = fj{l,q) = fj{li,l2,h,qi,q2,q3), i = i,2,3. 

Because the Ij are independent and in involution, it follows that dfj/dpk = dfk/dqj 
on the level surfaces {Ij = Cj, j = 1, 2, 3}. Therefore there is a generating function 
S{q, I) such that 

(7.9) dS/dqj = fj{I,q)=Pj, J = 1,2,3. 
Putting Qj = dS/dlj, it follows that 

dS = T^pjdqj + 'EQjdIj 
and therefore 

(7.10) in = J^dpj A dqj = J^dlj A dOj. 

By (7.5) the functions Ij are real on SU{3) and we may replace the Gj in (7.10) by 
their real parts (if necessary) to obtain real Darboux coordinates. In the remainder 
of this section we show that the angle variables Gj are elliptic functions. 

We turn to equation (7.8). We already have pj = loglj, j = 1,3. The re- 
maining equation can be obtained, in principle, once we have a nontrivial identity 



27 

Proposition 7.2. Let ( = e^^ = aiia22tt33, -^o = -^2 — -^ i — -^3 + 1- The following 
identity holds on 5'L(3, C); 

(7.11) (C/o + 2/1/3)2 + 4 cos2(1q2)(C - /i)(C - /3)(C - hh) = 0. 

Proof. We make extensive use of the identities defining the cofactors: A12 = 
^23031 — ^21^33 and so on, and of the corresponding identities coming from the 
inverse matrix A : ai2 = A23A31 — A21A33 and so on. In the following computation 
each term in braces is replaced by its expression obtained from such identities in 
order to pass to the next in the sequence of identities. 

I2 = 022^22 = ^11^22^33 — {fl22fll3}{fl22ft3l}/fl22 

= aiia22a33 - [{^13^31} + {ai2a2i}{a23a32} - (012^13^23 + a2iA3ia32)]/a22 
= /i + /s + 1 - 2AiiAis/a22 + (012^13^23 + a2i^3ia32)/a22 

or 

(7.12) ai2Ai3a23 + 021^31^32 = a22/o + 2A11A33 

= (C/o + 2/1/3) /an 033. 

On the other hand 

(7.13) 

(ai2Ai3a23)(a2i^3ia32) = (ai2a2i)(^i3^3i)(a23a32) 

= («iia22 - ^33)(^ii^33 - a22)(a22a33 - ^11) 

= — (C-/3)^— (/1/3-C)— (C-/i) 

a33 «iia33 «ii 

Let / = 012^13023 and g = 021^31032. Then 
(7.14) 

(/ + gf = fgiv^g + ^Wf? = fg[^M-\m) + eM\i<i2)? 

= 4fgcos'^{-q2). 

Combining (7.12), (7.13), and (7.14), we obtain (7.11). 
Returning to the generating function 5", we have 

dS/dqi =pi= log /i, dS/dqs = ps = log /s, dS/dq2 = P2 = log C 



so 



S(I, q) = Qi log /i + 03 log h + / log C(/, u)du. 



f 



q2{C,i) 
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The angle variables corresponding to the action variables Ij are 

ds n^ 1 dC 



dij I J J cdij 



dh J C dlf""- 
We rewrite (7.11) in the form 

(7.16) $(C, /, u) = F(C, /) + cos2(V2)G'(C, /) = 

to define ^ or w in terms of tt or ^ and /. Thus 

dC __d^/dlj _ {d^/dlj){d^/du) 
dfj ~ ~ d^/dC ~ dC/du ■ 

In particular, on the surface $ = 0, 

cos(^«) =i[F/G]^; sin(lti) = [1 + F/G]^; 
so 

ds n^ 1 ac , /"^^^'^"^ 1 

This can be written as a Jacobi elliptic integral; cf. [Co, p. 400]. Set z^ = 
{( — a)/{( — (3) where ct, (3 are the nonzero roots of F + G; then the last integral 
becomes 

(7.17) 62 = ^ / ^=, A;2 = /?/«. 

Note that a/3 = /1/2/3 is real and positive on 5't/(3). 

The other angle variables are also elliptic integrals. Straightforward calculation 

gives 

The first of the three integrals in each line are equal to —©2 and the remaining 
two transform, under the same change of variables as in (7.17), into sums of Jacobi 
elliptic integrals of the first and third kinds. 

To complete our discussion of SU{3) we consider the integration of the fiow with 
Hamiltonian I2, in the original coordinate system. The Poisson bracket determined 
by the foliation of the form zO differs from the Poisson bracket (5.11) by a factor 
—i. Thus the fiow on SU{3) is given by 

(7.18) / = -2(log/2,/) = -z(log 022^22,/), / e C^iSLin)). 
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Theorem 7.3. On SL{3,C), let 

p = (11/2/3)2, iV = — log(aiia22a33Mll^22^33)- 

2i 
These functions and the functions ajkAjk and ^log{ajk/Ajk) are real on SU{3). 
Under the flow (7. 17) uj evolves according to the pendulum equation 

(7.19) to = —2psmuj. 

Moreover 

(7.20) each ajkAjk is an algebraic function of cos uj and the Ij; 

(7.21) each time derivative of [^log{ajk/Ajk)] is an algebraic function of cos u 
and the Ij. 

Proof. A direct but somewhat tedious calculation using (7.18), (5.11), and various 
identities for cofactors gives 

(i) = — log(a22M22) "■ = i{aua22a33 - ^11^22^33) 
= — 2psinc<;. 

To obtain (7.20) we use the identities 

ttjiAki + aj2Ak2 + ttjsAks = Sjk = aijAik + a2jA2k + asjAsk 
which come from (7.05) to obtain 

djkAjk + ttkjAkj = 1 — Ij — Ik + Iij 
for distinct j, /c, £. Also 

djkAjkttkjAkj = ajkttkjAjkAkj = (ajjUkk — A£){AjjAkk — au) 
= Ijik + Ii — 2p cos UJ. 

Therefore ajkAjk and akjAkj are the roots of a quadratic equation with coefficients 
which are polynomials in the Ij, p = (/1/2/3) 2 , and cosuj. 

Finally, we consider log(ai2/Ai2). Another direct calculation gives 

= -j{aiia22{AiiA22 - 033) + ^11^22(011022 - ^33)} 

_ 1 p A11A22AS3 + 011022033 - 2/1/2 -, 
~ 4^ ai2Ai2 

■(pcoscu - /1/2). 



2012A21 

The calculations for other such terms are similar, and (7.21) follows, using (7.20). 
Note that on SU{3), the functions considered above are ajkAjk = |ojfcP and 
A lofffa-ii-Moi.) = are a^i.. 
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8. General nondegenerate J 

In this section we discuss the case of a spectral problem (3.2) whose characteristic 
matrix J has n distinct eigenvalues but is not otherwise constrained. Thus we may 
assume 

(8.1) J = diag(zAi, zA2, . . . , «An), Aj's distinct. 

The corresponding space of potentials P, the 2-form Op, and the associated Poisson 
bracket { , }p are defined as in section 2. The (continuous) scattering data s or 
(t;_i_, V-) which correspond to a potential g in P is a matrix- valued function or pair 
of such functions, defined on the set 

(8.2) T. = {^eC:Re{i\j^)=Re{i\kO, some jV A;}. 

This set is a union of lines through the origin; we will consider it as a union of 
rays from the origin and orient each ray from to oo. We will describe the spaces 
SD = {s} and SD' = {{v-^.,V-)} in more detail below. The analogue of Theorem 
3.1 carries over to this more general setting; [BCl]. Therefore the 2-form and 
Poisson bracket can be carried over to a form O5 and Poisson bracket { , }s on 
scattering data. In this section we prove the existence of Darboux coordinates and 
the complete integrability of the linear flows. 

Theorem 8.1. There are functions p^,qi,,l < z/ < (n^ — n)/2, holomorphic on 
a dense open algebraic subset of SL{n), which have the following properties. Let 
El, ... , E^ be the rays of E. There is an assignment of rays v 1-^ E^^^) such that 

(8.3) s in SD is uniquely determined by the values 

K(s(0, qMO). ^ e Efc(,), 1 < z/ < (n^ - n)/2}. 

(8.4) 05 = E/s,,,/P.A5g.. 

(8.5) {p^,c,P^,r;}s = {qfA,c,qu,v}s = 0; {p,.,^,qu,v}s = S^J{C - r]), ^,1] G 

(8.6) For any f in C°°{SL{n)), the distributions {f^,Pu,r]} have support at ^ = 
V e Efc(^), all ly; as before, f^{s) = f{s{^)). 

Theorem 8.2. The functions p^, of Theorem 8.1 may be chosen so that for each 
traceless diagonal matrix n, the Hamiltonian for the flows (2.8), (3.20) is a linear 
combination of the functional 



f ^"pAsmdt 
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The machinery needed for the proofs of Theorems 8.1 and 8.2 has aheady been 
developed in sections 3 and 4. To show how it apphes, we must describe the space 
of scattering data. Assume / [[^(x)!!^^ < 1 and consider the spectral problem (3.1), 
(3.2). Again there is a unique solution •0(-,2). This solution is holomorphic with 
respect to z^z G C\S, and its boundary values satisfy 

(8.7) ^{x,{\^%^)i) = ^{x,{:\-%^)i)v{i), ees\o. 

Given ^ in S\0, let H^ : M^ — > M^ be the projection defined by 

ttj-fc if Re{i\ji) = Re{iXkO, 



Then the limit 



{Ilfa)jk 

0, otherwise. 



(8.8) hm U^i;{x,Oe-^^' = siO. 

a;— >+cxD 

exists on E\0. There are factorizations 

(8.9) viO = v.{0-W(0 = s.iO-'s+iO; «±(0 = s{Ov±{0- 
The factors in (8.9) are characterized by the following conditions: 
(8.10) 

n5t'±(0 = ^±(0; n^«±(6 = «±(0; diag v± = i; 

if Re i{Xj — Xk)w > for w = {1 + ie)^, small £ > 0. 

There are further conditions on s and v as functions of ^ G E, which we do not 
need to cite here; see [BCl] for a complete discussion of these conditions and of the 
algebraic facts we are using in this section. 

Proposition 3.2 carries over, in the following form [BC3]: 



^.11) ns = — [ tr(v-(5v)vT^ A s-^5s)dz. 

2m J-^ 



The strategy for proving Theorems 8.1 and 8.2 is the same as for proving Theo- 
rems 3.7 and 3.8: pointwise analysis of the form O5. 

Suppose ^ is in E\0. After conjugation by a permutation matrix (which depends 
on the ray of E containing ^), we may assume that 

(8.12) Re iXiC = ■ ■ ■ = Re iXd,^ > Re iXd,+iC = ■■■ 
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Then v{^) and s(^) have the block diagonal form 



^.13) 



Ai 
A2 








A. 



Aj e SL{dj 



Moreover V-{^) and s+(0 have this form and are upper triangular, while t'+(^) and 
S-{^) have this form and are lower triangular. With this normalization, the point- 
wise 2-form being integrated over the ray containing ^ is (a constant multiple of) 
the sum of the forms as in Theorem 4.1 for the matrix groups GL(di), GL{d2), .... 
There is an analogous decomposition of the Poisson bracket, which is computed as 
in section 3. 

It follows from these considerations that the (trivial) extension of the results in 
sections 4 and 5 to block diagonal matrix groups yields the desired functions pi,, q^ 
of Theorems 8.1 and 8.2. 
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